Rotational Doppler effect in left-handed materials 
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We explain the rotational Doppler effect associated with light beams carrying with orbital angular 
momentum in left-handed materials (LHMs). We demonstrate that the rotational Doppler effect 
in LHMs is unreversed, which is significantly different from the linear Doppler effect. The physics 
underlying this intriguing effect is the combined contributions of negative phase velocity and in- 
verse screw of wave-front. In the normal dispersion region, the rotational Doppler effect induces a 
upstream energy flow but a downstream momentum flow. In the anomalous dispersion region, how- 
ever, the rotational Doppler effect produces a downstream energy flow but a upstream momentum 
flow. We theoretically predict that the rotational Doppler effect can induce a transfer of angular 
momentum of the LHM to orbital angular momentum of the beam. 

PACS numbers: 42.25.-p; 42.79.-e; 41.20.Jb; 78.20.Ci 
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I. INTRODUCTION 



After the first experimental observation of negative re- 
fraction [l|, Q , some intriguing or even counter-intuitive 
phenomena in left-handed materials (LHMs), such as am- 
plification of evanescent waves [1, 0] , unusuaLphoton tun- 
neling [1, 0] , negative Goos-Hanchen shift [J [1] , and re- 
versed linear Doppler effect [§, have attracted much 
attention. Linear Doppler effect is a well-known phe- 
nomenon by which the frequency of a wave is shifted 
accordin g to the relative velocity of the source and the 
observer [ll|. The conventional understanding of linear 
Doppler effect is that increased frequencies arc measured 
when a source and an observer approach each other. Ap- 
plications of the effect are widely established and include 
radar, laser cooling, flow measurement, and the search for 
astronomical objects. The inverse linear Doppler effect 
refers to frequency shifts that are in the opposite sense 
to those described above. For example, increased fre- 
quencies would be measured on reflection of waves from 
a receding boundary d, [l3| . 

The interaction of a moving medium with a Lagurre- 
Gaussian (LG) field raises the fundamental question of 
the Doppler effect. As a medium moves across the heli- 
cal wave fronts of the LG field, it experiences, in addition 
to the usual linear Doppler shift related to the veloc- 
ity in the light propagation direction, a most intriguing 
frequency shift, the so-called rotational Doppler effect 
associated with the azimuthal velocity. The rotational 
Doppler effect results in frequency shift when the light 
beam is rotated around its propagation axis. The rota- 
tional Doppler effect has been extensively studied in reg- 
ular right-handed materials (RHMs) [H, [H Hi, [H, [11. 
Veselago has suggested that LHMs would reverse nearly 
all known optical phenomena ^17j . Now an interesting 
question arises: what happens to the rotational Doppler 
effect in LHMs? When the LG beams incident into 



LHMs, the helical wave fronts should reverse their screw- 
ing fashion [l^ . Thus we can predict that the rotational 
Doppler effect will result in many intriguing phenom- 
ena. The investigation of the rotational Doppler effect 
will provide insights into the fundamental properties of 
LHMs and will allow us to better understand the inter- 
action of light with LHMs. 

In this work, we try to reveal the rotational Doppler ef- 
fect associated with light beams carrying orbital angular 
momentum in LHMs. First, starting from the Maxwell's 
equations in moving frame, we obtain the analytical de- 
scription for LG beams propagating in rotating LHMs. 
Our formalism permits us to introduce the effective in- 
dex to describe the wave propagation. Next, we attempt 
to recover how the wave-front evolves, and how the screw- 
ing wave-fronts result in upshifted and downshifted fre- 
quencies. In order to explore the rotational Doppler ef- 
fect, we examine a simplest intensity pattern formed by 
LG beams with equal magnitudes and waist parameters. 
Then, we want to investigate how the rotational Doppler 
effect gives rise to the anomalous upstream transverse 
Poynting vector and angular momentum flow. Finally, 
we attempt to explore how the rotational Doppler ef- 
fect influences the rotation of intensity pattern inside the 
LHM. 



II. EFFECTIVE REFRACTIVE INDEX 

To investigate the effective index of rotating LHM, we 
firstly use the Maxwell's equations to determine the field 
distribution inside the LHM. We consider a monochro- 
matic electromagnetic field, E(r, t) ~ Re[E(r) exp{—iujt)] 
and B(r,t) = Rc[B(r) exp(— icji)], of angular frequency 
uj propagating from a conventional RHM to the LHM. 
The field can be described by Maxwell's equations 

dB 

VxE = - — , V-B = 0, 
ot 
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V X H = 



dB 



V • D = 0. 



(1) 
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From the Maxwell's equations, we can easily find that the 
wave propagation is permitted in materials with e < 
and /i < 0. Veselago termed these LHMs, because the 
vectors E, H and k form a left-handed triplet instead of 
a right-handed triplet, as is the case in regular RHMs 

In order to investigate the effective index, we introduce 
a rest frame and a moving frame fl^. In the rest frame 
the LHM rotates with an angular velocity v,^ — ft x r 
and in the moving frame the LHM is at rest. We restrict 
our analysis to a small velocity with v c. For the LHM 
at rest, we have the following constitutive relations: 



Q> 



D' = eoe{uj')F,' B' = ^oAi(w')H' 



(2) 



where we have used primes to denote the fields and their 
frequencies in the moving frame. The fields in the moving 
frame can be expressed in the rest frame by a Lorentz 
transformation, which gives to first order in v/c: 



D 



B 



V X H 

V X E 



ene(t-')(E + v X B), (3) 
AioMM(H- V X D). (4) 



Here electric permittivity e and magnetic permeability fi 
are given as a function of the frequency in the moving 
frame. We also assume that e and /i depend only on the 
frequency and is independent of the state of motion of the 
LHM. On combining these two equations, we can obtain 
the following wave equation: 

V^D ^ ^n{Lo'^fuj^D + 2[n{uj'^f~l]nLu{a±d^)D, (5) 

where we have introduced the relation n{oj') ~ 
— £{ui')^,{uo'), and a — ±1 corresponds to left- and 
right-handed circularly polarized light, respectively. To 
determine the effective indices, we solve the wave equa- 
tion for helical wave- front LG beams. 

We introduce the Lorentz-gauge vector potential to de- 
scribe the propagation characteristics of LG filed. The 
vector potential of the beam propagating in the +z di- 
rection can be written in the form 

A = j4o(aej, + f3ey)up,i(r) exp(mfcoz — icut). (6) 

Here Aq is a complex amplitude, n is the effective index, 
ko = u>/c, c is the speed of light in vacuum, and ey 
are unit vectors, respectively. The coefficients a and /3 
satisfy the relation a = i{a(3* — a*/3). It is well known 
that each photon in a light beam carries a spin angular 
momentum. 

To be uniform throughout the following analysis, we in- 
troduce different coordinate transformations z*(i — 1,2) 
in the RHM and the LHM, respectively. The field can be 
written as 



pi 



X exp 



X exp[ 



V2r 
—ink^r^z* 



exp[±iZ(^] 



exp 



w^{zl) 



-i(2p + ;-f l)arctan(z*/z_R)], (7) 




(b) 



FIG. 1: (Color online) The helical wave front for Laguerre- 
Gaussian beam result from an azimuthal phase structure of 
exp[±ii</5]. The LHM is rotating with angular velocity Q,. 
The phase velocity Vp reverses its direction and the wave- 
fronts reverse their screwing rotation in the LHM. When a 
metamaterial "atom" moves across the helical wave fronts, it 
experiences the so-called rotational Doppler effect, (a) The 
helical wave front with exp[-|-i<p] and the effective frequency 
uj' = u}-\-Q,. (b) The helical wave front with exp[— and the 
effective frequency w' = w — n. 



w{zl) = w^^1 + {z*Izr)\ R{z*) = z* + 



(8) 



Here Cpi is the normalization constant, Lp[2r^/w^(z*)] is 
a generalized Laguerre polynomial, zr = nfcoWo/2 is the 
Rayleigh length, and w{z*) is the beam size and R{z*) 
the radius of curvature of the wave front. The last term 
in Eq. ((7|) denotes the Gouy phase which is given by 
$ = -(2j3 + / + 1) arctan(z*/zfl). 

Now let us study the screw of the wave front. Here the 
sense of the positive angles is chosen as anticlockwise, 
while negative angles are considered in the clockwise di- 
rection. In the rotating LHM, the constant wavefront 
satisfies 



n+koz* + ■ 
n-koz* + 



-in-kQT^ 



+ lip + *+(z*) = const, (9) 
-^V? + const. (10) 



The schematic view of the wave front is a three- 
dimensional screw surface of (rcos(p, rsinip, z). The 
plotting range of r is from to 5wo with the interval of 
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Ar = 0.5wo and that of |n±|fco2; is from to iw with 
the interval of |n±|A:oAz = O.Itt. The wavefront struc- 
ture reverse their screw types with a pitch of Ao/|n+ 1 and 
Ao/|n_|, respectively (see Fig.[T]). 

The helical wave fronts with clockwise or anticlockwise 
screwing fashion will result in different effective indices. 
Substituting the LG beam in the wave equation, we have 

4M - n^c.'^)"^ - 2[n'{c.'^) - ± 0- (H) 

Thus from Eq. pip , and using the relation k(uj) = 
n(u!')uj/c, we can obtain an equation of the refractive 
indices associated with the waves of each screwing fash- 
ion: 



n±{uj) = n{uj'j_) 



1 



-{<j±l). (12) 

UJ 



Here we have approximated the refractive indices to the 
first order in fl/uj. To obtain the indices, we need to know 
the effective frequency. Strictly speaking, we cannot ob- 
tain it from the wave equation in the moving frame. Now 
a question arises: how to determine the value of effective 
frequency in the LHM? 



III. ROTATIONAL DOPPLER EFFECT 

In the case of electromagnetic radiation this usually 
means that the subunits must be much smaller than the 
wavelength of radiation. Then the unit cells of metama- 
terials can be modeled as the atoms (or molecules) in 
ordinary materials. In particular, electric atoms with an 
electric-dipole moment leading to a negative electric per- 
mittivity, and magnetic atoms with a magnetic-dipole 
moment leading to a negative magnetic permeability. 
Hence the rotational Doppler effect in LHMs is the result 
of the optical properties of the individual metamaterial 
"atoms" (see Fig. [T]). 

To quantify the changes that occur in the Doppler 
effect for metamaterial "atoms" interacting with light 
beams with orbital angular momentum, the quantum me- 
chanical derivation [l3| is developed for describing the 
rotational Doppler effect in LHMs. The Doppler shift 
experienced by a LG field in a moving frame is given by 



{2p + l + l)zii 

^! + 4 



+ 4) 



2z? 



nkor 

1 

R 



1 



±1 



where Vz, Vr, and are the axial, radial, and azimuthal 
velocity components of the metamaterial "atom" , respec- 
tively. The axial Doppler shift is obtained from nk^Vz- 
There are two additional terms which account for the 
radial phase and the Gouy-phase. Note that the two ad- 
ditional terms should be reversed, since the negative in- 
dex and the inverse Gouy-phase shift [l^]. Thus we can 
conclude that the linear Doppler effect associated with 



light beams is slightly different form the counterpart of 
plane wave |1_0,] . Further investigation shows that the 
radial Doppler shift, given by the term proportional to 
Vr, should also be reversed in LHMs. 

A most intriguing frequency shift, the so-called rota- 
tional Doppler effect arises from the azimuthal velocity 
= U,r. The rotational Doppler effect results in fre- 
quency shift when the light beam is rotated around its 
propagation axis. The corresponding rotational Doppler 
shift in LHMs is given by 5^q = (cr ± VjVl, which is pro- 
portional to the total angular momentum. It can be seen 
that beams propagate differently in a rotating LHM, de- 
pending on whether the wave fronts turn in the same 
rotation sense as the LHM or in the opposite sense. The 
frequency of an incident LG beam whose wave-front is 
screwing as the LHM is downshifted to w' = cj— {a + l)^, 
whereas a LG beam of the opposite sense is upshifted to 
Lo' = uj — (cr—l)fl. Unexpectedly, the rotational frequency 
shift in LHMs is unreversed. 

When a LG beam incident into a LHM, the helical 
wave fronts reverse their screwing fashion [l^|. Why is 
the rotational frequency shift unreversed? By closely ex- 
amining the evolution of wave- front, we find the phase 
velocity also inverses its direction (see Fig.[T]). Hence the 
dynamic evolution between the metamaterial "atom" and 
wave-front remains unchanged. We conclude that the 
physics underlying the unreversed effect is collective con- 
tributions of negative phase velocity and inverse screw of 
wave- front. 

Though the rotational Doppler effect is unreversed, 
much more counter-intuitive phenomena will be caused 
in the LHMs. In order to explore the rotational Doppler 
effect, we will examine a simplest intensity pattern which 
is formed by LGo,+i and LGq.-i beams with equal mag- 
nitudes and waist parameters. In principle, any arbitrary 
light beam can be described by a superposition of LG 
modes '21'] . In the case of equal frequencies, such a com- 
bination is well known to be equivalent to an ordinary 
Hermite- Gaussian (HG) beam with a double-spot inten- 
sity distribution. We will explore the energy flow and 
angular momentum flow caused by rotational Doppler 
effect. 

First let us examine the propagation characteristics 
of energy flow, which is usually discussed by use of the 
Poynting vector. There has been considerable interest in 
orbital angular momentum of LG beams relating to az- 
imuthal component of Poynting vector [13, IH, [3] ■ The 
time average Poynting vector, S, can be written as 



S = iRe[(E^ 



E_) X (H++H_)*]. 



(14) 



We denote the radial, azimuthal, and axial components 
of the vector S in the cylindrical coordinates by the no- 
tation Sr, S^p, and Sz, respectively. The component Sr, 
relates to the spread of the beam as it propagates. The 
azimuthal component S^p describes the energy flow that 
circulates around the propagating axis. The axial compo- 
nent Sz describes the energy flow that propagates along 
the -\-z axis. 
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FIG. 2: (Color online) In the normal dispersion region , the 
rotational Doppler effect induces (a) the downstream energy 
flow and (b) the anomalous upstream transverse momentum 
flow. The parameters are u = 1.366a;o and = 0. The green 
arrows denotes the transfers Poynting vector and momentum 
flow. Please note that the beam divergence is not shown due 
to normalization of transverse dimension. 



Next we attempt to describe the evolution of momen- 
tum flow. The momentum flow for the electromagnetic 
wave can be written in the form 11 ll 



T = -Re[(D • E* 
2 



-B-H*)/- (DE* +BH*)]. (15) 



Here T is the momentum flow also referred as the 
Maxwell stress tensor, and / is 3 x 3 identity matrix. 
For a HG beam, E = E+ + E_ and H = H+ + H_. 
Note that the momentum flow is signiflcantly different 
from that obtained by Minkowski or Abraham momen- 
tum [ll'|. The cross product of this momentum density 
with the radius vector yields an orbital angular momen- 
tum flow. The orbital angular momentum flow in the z 
direction depends upon the component of , such that 
= r X T^. 

In order to accurately describe the Poynting vector and 
momentum flow, it is necessary to include material dis- 
persion and losses [1^. Thus, a certain dispersion re- 
lation, such as the Lorentz medium model, should be 



introduced. The constitutive parameters are 



mp 



(16) 
(17) 



To avoid the trouble involving in a certain value of fre- 
quency, we assume the material parameters are ujeo = 



mp 



2u}^ and 7e 



7" 



0.25cjo. It 



should be noted that, as we will see in the following, the 
dispersion and negative indices will play a very important 
role in the rotation of Poynting vector and momentum 
flow. 

Figure [5] shows the characteristics of Poynting vector 
and momentum flow. The circumstance is really intrigu- 
ing: the novel azimuthal energy flow and momentum flow 
present in the HG fleld. Counter-intuitively, the angular 
velocity fl of the rotating LHM and the circulation of 
Poynting vector are opposite, as shown in Fig.[2Ka). The 
presence of the azimuthal momentum flow has been 
identified as the orbital angular momentum. We thus 
predict that the rotational Doppler effect can induce a 
transfer of angular momentum of the LHM to obitual 
angular momentum of the beam. The sign of the orbital 
angular momentum coincides with handedness of rota- 
tion of the LHM. Positive orbital angular momentum is 
associated with the counter-clockwise rotation > 0, 
as plotted in Fig. (Hb). On the contrary, the sign of 
the orbital angular momentum in the HG field is nega- 
tive: when the LHM rotates clockwise 17 < 0, the angular 
momentum flow seemingly corresponds to the clockwise 
circulation. The physical origin of transfer of angular 
momentum is rooted in the geometric phase difference 
between LGo,+i and LGo,-i beams. 

Whether the rotational Doppler effect induce the 
downstream or upstream angular momentum depends 
on the normal or anomalous dispersion. In the nor- 
mal dispersion region, dKe[n{uj)] / duj > 0, the rota- 
tional Doppler effect induces a upstream energy flow 
[see Fig. [2Ia)] but a downstream momentum flow 
[see Fig. [IJb)]. In the anomalous dispersion region, 
dRe[n{uj)]/duj < 0, however, the rotational Doppler ef- 
fect induces a downstream energy flow and a upstream 
momentum flow. Interestingly, we find that the trans- 
verse energy flow in the LHM is antiparallel to the trans- 
verse momentum flow. In our opinion, the main reason 
for such an inconsistence is the negative refractive index. 
Note that the novel phenomenon is significantly different 
from the counterpart in regular RHMs (26i |. 

To understand a better physical picture, let us reveal 
the instant field distribution in the LHM. The rotation of 
the intensity pattern in LHM arises from the difference 
in the refractive indices for constituent LG beams. The 
transverse beam pattern is determined by the geometric 
phase difference: 



Akz* + 



R+iz*) R.{z*) 



2lip - A$ = const, (18) 
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FIG. 3: (Color online) Instant spatial structure of HG beam 
screws with pitch Zp — 47r/(fc+ — fc_) in the LHM. The rota- 
tional Doppler effect induces rotation of intensity pattern: (a) 
the downstream rotation {il > 0, vj > 0) in normal dispersion 
region, (b) the upstream rotation (il > 0, < 0) in anoma- 
lous dispersion region. Please note that the beam divergence 
is not shown due to normalization of transverse dimension. 

where Ak = — k^, k± — n±{uj')uj / and A<i> — 
(f>+ — (f>_. The schematic view of the spatial structure 
of instant field is plotted in Fig. [31 Here we have ne- 
glected the phase difference caused by radial and Gouy 
phases. The spatial structure exhibits an screwing type 
with a pitch of Zp = 47r/(fc+ — k-) along the +z axis. 

In this case, the constituent LG beams propagate dif- 
ferently in the LHM which leads to a change in their 
relative phase. The pattern rotates with an angular ve- 
locity n7 — c{k+ — fc-)/2. The angle per unit length by 
which the intensity pattern is rotated is called the spe- 
cific rotary power A6'i = 27r(n_|_ — n_)/A. Hence the 
rotational Doppler effect induces a downstream or a up- 
stream rotation of the intensity pattern, depending on 
whether the frequency locates in the region of normal or 
anomalous dispersion. In the normal dispersion region, 
A;+ > the angular velocity > 0, and the rotational 
Doppler effect induces a downstream rotation of inten- 
sity pattern [see Fig. [3]Ja)]. In the anomalous dispersion 
region, k+ < k+, the angular velocity < 0, and the ro- 
tational Doppler effect induces a downstream rotation of 
intensity pattern [see Fig. [n^b)]. A further point should 



be mentioned that the large rotation can be obtained if 
and fc_ have opposite signs (i.e., fc+ > 0, fc_ < or 
A;+ < 0, fc_ > ). Our result suggests that the inter- 
esting images rotation [23, [H, [2^ might more readily be 
accessed in LHMs. 

The possibility of having an upstream rotation of in- 
tensity pattern in the presence of anomalous dispersion, 
i.e., in the presence of negative group velocity. As nega- 
tive group velocities in dispersive media are forbidde n by 
Kramers-Kronig relations in all frequency regions □J. 
Thus negative group velocities are possible onl y in the 
presence of absorption. In a recent experiment [30|, the 
negative group velocity have been demonstrated by ob- 
serving that the pulse advances in time with respect to 
the same wave packet propagating in vacuum. The low- 
loss LHMs are therefore good candidates for the exper- 
imental observation of significant upstream rotation of 
intensity pattern. Conversely, the anomalous upstream 
rotation of intensity pattern could provide an interesting 
way to reveal the nature of LHMs. 



IV. CONCLUSIONS 

In conclusion, we have investigated the rotational 
Doppler effect associated with light beams carrying or- 
bital angular momentum in LHMs. Unlike the reversed 
linear Doppler effect, the rotational Doppler effect in 
LHMs is unreversed. We have revealed that the physics 
underlying this intriguing effect is combined contribu- 
tions of negative phase velocity and inverse screw of wave- 
front. Though the rotational Doppler effect is indepen- 
dent of the refractive index, much more counter-intuitive 
phenomena will be caused in the LHMs. In the normal 
dispersion region, the rotational Doppler effect induces 
a downstream energy flow but a upstream momentum 
flow. In the anomalous dispersion region, however, the 
rotational Doppler effect induces a upstream energy flow 
but a downstream momentum flow. We theoretically pre- 
dict that the rotational Doppler effect can induce a trans- 
fer of angular momentum of the LHM to orbital angular 
momentum of the beam. It is possible that the study 
of Poynting vector and momentum flow in LHMs may 
make a useful contribution to long established Abraham- 
Minkowski dilemma. We have shown that the momentum 
flow in LHMs should be antiparallel to the Poynting vec- 
tor. Thus the direction of momentum flow is signiflcantly 
different from that obtained by Minkowski or Abraham 
momentum. We consider that this difference might help 
us to illustrate this long established problem. 
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